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We study theoretically the stability of flow in superfluid 3 He-A. The calculations are done using a 
one-dimensional model where the order parameter depends only on the coordinate in the direction 
of the superfluid velocity v s . We concentrate on the case that the external magnetic field H is 
perpendicular to v s , where only few results are available analytically. We calculate the critical 
velocity v c at which the superflow becomes unstable against the formation of continuous vortices. 
The detailed dependence of v c on the temperature and on the form of the underlying orbital texture 
l(r) is investigated. Both uniform and helical textures of 1 and two types of domain-wall structures 
are studied. The results are partially in agreement with experiments made in a rotating cylinder. 

PACS: 67.57.Pq 



I. INTRODUCTION 

The superflow of 3 He-A differs markedly from the well 
known superfluid 4 He, where the decay of persistent cur- 
rents is prevented by topological constraints. The circu- 
lation of the superfluid velocity v s on a closed contour is 
not quantized in 3 He-A, but depends on the field l(r) of 
the orbital anisotropy vector. Thus 3 He-A can respond to 
externally applied flow by forming an inhomogeneous tex- 
ture of l(r). The rigidity of the order parameter stabilizes 
the uniform bulk state (with l(r) = constant) for small 
flow velocities in most cases, depending on the magnitude 
and direction of the external magnetic field H. With 
increasing velocity this configuration becomes unstable 
against textural inhomogeneities that finally lead to the 
formation of "continuous" vortices. In several cases a 
helical texture of l(r) is stabilized at intermediate veloc- 
ities. 

The flow properties of bulk 3 He- A were studied in 
many theoretical papers which peak in the late 1970's 

p2fl . The majority of that work studied the case where 
the superfluid velocity and the magnetic field are parallel. 
Various methods were used to generate the flow experi- 
mentally (torsional oscillator, thermal gradient, a piston, 
etc. see Refs. fl3| ) fl4| for reviews). Although some of the 
predicted features were seen in the experiments, no sat- 
isfactory agreement between theory and experiment was 
found. 

More recently, a rotating cryostat has been used to gen- 
erate superflow . This method gives a well defined dc 
superfluid velocity under very steady conditions. These 
experiments motivated our theoretical studies. Firstly, 
it was necessary to study the case where v s is perpen- 
dicular to H. The helical texture in this case was pre- 
viously considered only in one paper 0]. Secondly, the 
calculation has to be generalized to temperatures sub- 
stantially below the superfluid transition temperature T c . 
Thirdly, we consider two different types of "solitons". 
These are domain-wall like structures whose effect on the 



flow properties were previously considered in a few cases 
P,pUl3-|l3j . Here we present detailed calculations of the 
critical velocities v c for the appearance of helical textures 
and of vortices in different cases with v s _L H. The com- 
parison of the results with experimental measurements 
has been given before fill . 

The hydrodynamic free energy describing the current- 
carrying states of 3 He-A is recalled in Sec. H In Sec. 
Ill we study the general features of our one-dimensional 
model and its numerical solution. Sec. |v| discusses the 
instabilities of uniform and helical textures. In Sec. |v] 
the critical velocity associated with two initially inhomo- 
geneous textures, a dipole-locked and a dipole-unlocked 
soliton, is investigated. The comparison to experiments 
is briefly discussed in Sec. VI. 



II. HYDROSTATIC THEORY 

The order parameter of 3 Hc-A is fully specified by 
defining an orthonormal triad {rh, n, 1} and a unit vector 
d. The triad describes the orbital part and d the spin 
part of the tensor order parameter p6| 



A M = Adn(mj + ihj). 



(1) 



Here A is a (temperature-dependent) constant and 1 = 
m x n. All the unit vectors may vary as functions of 
location r. A change in the total phase <E> of the order 
parameter ([!]) is given by 

5 * = \ [™i 6fl 3 ~ ^jSrhj] = rhjShj. (2) 



This allows for the definition for the superfluid velocity 
as 



2m 



(3) 
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where m in the prefactor is the mass of a 3 He atom. The 
equilibrium properties of the superfLuid are determined 
by the order-parameter configuration that minimizes the 
total free energy. In the hydrodynamic approximation it 
has the form 



F = / d 3 r / = 



/ d 3 r (/ gr 



+ /d + /h). 



(4) 



The gradient energy density can be written as 

/gr = \P^I + f - P±)(l ■ V s) 2 

+ cv s -vxi-c (i-v s )(i-vxi) 

+ \K S {V • i) 2 + \K t (\ • V x i) 2 + \K h \\ x (V x i)| 2 
+ \Ks\(\ ■ V)d| 2 + iif 6 Ey[(i x V)^,-] 2 . (5) 

The first four terms give the kinetic energy of the 
anisotropic supcrfluid. The C and Co terms give cou- 
pling between flow and an inhomogeneous 1 field. The 
remaining five terms are the bending energy densities for 
1 and d. 

The energy density of dipole-dipole interaction is 



/d 



(6) 



Comparing this to the kinetic energy (|J) defines the 
dipole velocity — */<?d/p|| ~ 1 mm/s and the dipole 
length £d = h/2mvd ~ 10 /im. In addition, the energy 
density in the presence of an external magnetic field H 



/ h = i<? h (d.H) s 



(7) 



It is customary to define the dipole field iJd = \J ffd/ffh ~ 
2 mT by comparing (JsJ) and (jjj). 

All the coefficients p±, P\\, C, Co, K s , K t , Kb, K 5 , K e , 
g&, and gh are positive and depend on the temperature 
T and the pressure p. Their values are determined as 
explained in Ref. |l7f . In particular, the coefficients are 
calculated using consistently the weak-coupling approx- 
imation. In reduced units of i>d, £d and our results 
are independent of ga and gh, but they are needed for 
comparison to experiments |lq] . For g d we write 



g d (T,p)=4g a D (p) A 2 A (T,p), 



(8) 



where Aa is the maximum energy gap in the weak- 
coupling approximation. All calculations are done at the 
melting pressure where g^ = 5.9 • 10 44 J _1 m -3 Mj. We 
wish to emphasize that the calculations contain no ad- 
justable parameters. 

The hydrostatic theory gives a good description of the 
superflow in 3 He-A over most of the temperature region 
< T < T c . It becomes invalid in a small region T c — T < 
10 T c around T c , where the Ginzburg-Landau critical 
velocity ~ ^T^gl ls smaller than v&. We also limit to 
such low fields that the deformation of the A phase order 
parameter (|l]) towards the Ai phase can be neglected. 



III. ONE-DIMENSIONAL CALCULATION 

We study flow in bulk 3 He-A far from any walls. The 
main assumption in the present work is that the order 
parameter ([!]) depends only on one spatial coordinate x. 
It follows from the definition (|^) that v s is always parallel 
to the x axis, v s || x. In addition to the homogeneous 
state, the ID model allows us to calculate the structure of 
helical textures and the deformation of solitons in a flow 
that is perpendicular to the soliton wall. Moreover, we 
also can determine the stability limits of such textures 
against ID perturbations. We will argue in section VI 
that the local stability of the states we consider is indeed 
determined by such perturbations. 

The cases of H = §10,011 and H || x | |,|§Jll 
have been studied extensively in the literature . Here we 
study the more complicated case H ± x [pH. We study 
in particular the high field limit H 3> H dl where dlH 
everywhere. 

The order parameter of 3 He-A can be parametrized by 
introducing three Euler angles a, (3 and 7 for the orbital 
triad {m, n, 1} and polar and azimuthal angles i/j and 
<f) for d. In this representation one has to beware the 
unphysical singularities at the poles sin [3 — or sin tp = 
0. For this reason we choose the polar axis z of the fixed 
coordinate system perpendicular to the direction of v s 
(and parallel to the external magnetic field). This results 
in a more complicated form of the energy functional but 
it enables us to avoid the singularities in all stationary 
states. With these definitions 

1 = sin (3 cos a x + sin (3 sin a y + cos (3 z (9) 
d = sin ip cos cj) x + sin ip sin (j> y + cos i/j z (10) 
ft ,dl „ da s „ 

The form of the unknown functions a(x), (3(x), 7(2;), 
tjj(x) and (j)(x) in equilibrium corresponds to the mini- 
mum of the free energy functional (Q) where 
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and similarly for other angles. 

We point out two analytic observations that are useful 
in testing the convergence and the accuracy of the nu- 
merical solution. Firstly, the superfluid velocity v s (p"l|), 
and therefore also the total free energy (^), do not de- 
pend explicitly on the angle 7 but only on its derivative. 
The corresponding Euler equation ([l2]) reduces to a con- 
servation law 



Of 



d{dj/dx] 



— —p = const. 



(13) 



The quantity (2m/h)p is the x component of the super- 
current density 
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P±_V S + 0|| - p±)l(l ■ v s 

C*V x 1 — C 1(1 - V x i). 



(14) 



Another conserved quantity in the problem is the one 
corresponding to the fact that the free energy ([|) does 
not depend explicitly on x either. The invariant related 
to this is analogous to the Hamiltonian of classical me- 
chanics and can be brought to the form 



,/gr - ./d - h = const. 



(15) 



We wish to study a one-dimensional interval of length 
L>(d. A flow through the system is achieved by keep- 
ing a fixed phase difference A$ = $(L/2) - L/2) 
between the endpoints of the line. The constancy of A<f> 
as a function of time t is enforced by imposing the bound- 
ary condition 



dA$ 



/ d7 ^,da\ 
= - -r- + cos/3— 
dt \dt ' dt J x= , 

= 0. 
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-r + cos/3— 

dt ' dtj x= _ 



(16) 



It is, however, more advantageous to express the re- 
sults in terms of a driving velocity defined by v n = 
(h/2m)A$/L. This is a more convenient quantity than 
A$ because all our results are independent of L when 
expressed in terms of v n . v n could be identified as the 
velocity of the normal component that drives the super- 
fluid component of the liquid. However, v n should be 
considered as a scalar parameter since all vector quanti- 
ties in this paper (like v s and j s ) are given in the frame 
where the normal fluid is at rest, v n = 0. (See Ref. ]l9| ] 
for a general formulation with v n ^ 0.) 

In general case we have to use numerical methods 
to determine the equilibrium state of the system. The 
order parameter is taken to be defined at N equally 
spaced discrete points on the line. The discretization 
length Ax = L/N is chosen much smaller than the dipole 
length, usually Ax < 0.1£d- As a first step we have to 
choose some initial configuration for the five angles. A 
given A<I> or v n is implemented by taking an initial guess 
j(x) — —^v n x. The angle functions are then iterated 
numerically towards the equilibrium state for a given v n 
using the following diffusion-like equations 

\x\ sin ft 
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together with the boundary condition (|l6|) . The simplest 
discretized expressions have been used in representing the 
derivatives in equations (|T7|). Since we do not attempt 
to describe the true time evolution of the textures, the 
viscosity constants /ii , /i2 and /Z3 can be chosen according 
to numerical convenience. 



IV. UNIFORM AND HELICAL TEXTURES 

The simplest texture has constant 1 || d || x. This 
uniform state minimizes both the dipole-dipole energy 
(0) and the field energy (0) for Hlx. It also corresponds 
to the minimum of the first two terms in the gradient 
energy (||) because p|| < p±. The current in this state is 
linear in v n , j s = p||W n x, as illustrated in Fig. |l|. 
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FIG. 1. The supercurrent j s<x as a function of driving 
velocity v n at T — 0.6T C and H S> The diagonal line 

corresponds to the current in the uniform texture, j s = p\\v n . 
The solid line between v c i and v C 2 corresponds to helical tex- 
ture at the optimal wave vector q op t- 



The uniform state is stable at small velocities (if H ^ 
0). We call the stability limit of the uniform texture 
as the first critical velocity v c i. At v n = v c i the uni- 
form state becomes unstable against a helical deforma- 
tion where 1 winds around the direction of the flow, see 
Fig. |^. At this point the energy cost in forming an in- 
homogeneous texture is compensated by reductions in 
other energy terms. In particular, the superfluid velocity 
v s ( fil| ) is lowered for a given v n and there is a negative 
contribution from the energy term with the coefficient Co 

(I) !• 

The instability point can be studied by expanding the 
energy (Q) around the uniform solution. This calculation 
was done by Lin-Liu et al in the Ginzburg-Landau region 
0. We generalize this calculation to all temperatures. 
We define 7 as above but otherwise use a parametrization 
that is different from (|9|-pT|): 
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FIG. 2. Variation of 1 in a helical texture in one spatial 
dimension (x is in the direction of flow). One wavelength A 
of the helix is shown. 



i = [i - Wv + - Wl + l D 2 \* + + w (is) 

d = [1 - \{dl + dl) - \{dl + dlf]± + d y y + d z z (19) 



s_2mi dx + 2[ly dx z dx )[ 



if/ 2 

4\ L y 



(20) 



These equations are valid up to fourth order in l y , l z , d y , 
and d z . We eliminate 7(2:) in favor of the current p (Jl^ ) 
by defining a new free energy G = F + J dx(dj /dx)p. 
Substitution of (|i"g|-p0|) into G and expansion to second 
order gives linear Euler-Lagrange equations. These have 
the solution 



ly{x) 

l z (x) 

dy(x) 

d z (x) 



us svn(qx) 
(u/s) cos(qx) 
uSiS sm(qx) 
(1162/ s) cos(qx) 



(21) 
(22) 
(23) 
(24) 



where q is the wave vector of the helix. The magnetic 
field in the transverse z direction introduces "easy" and 
"hard" directions for the amplitudes, and thus the helix 
has an elliptically distorted form. When G is minimized 
with respect to s, S±, and 62 we find 
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H 2 



61 = (K 5 q 2 
h = {K 5 q 2 
s 2 = K2/K, 
K? = l + (p x - l)p 2 + K h q 2 



For the free energy we find the expansion 
G = Gq 

where 
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B=[(2 P± -l)Cb-J] 
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FIG. 3. The stability of different states in the v n — q plane 
at T = 0.6T C and for H » H d . The solid line marks the 
instability of the uniform texture according to A(v n ,q) — 
(Jil|). The helical texture is stable inside the region spanned 
by the dashed line. The dash-dotted line shows the optimum 
wave vector q op t as a function of the driving velocity v n - 



+K 5 5l)q 2 3{p ± 1) V] 
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+(1 - (50(1 - 8 2 )(6i + 6 2 ) - 2{p± - iy P 



(32) 



For simplicity, we have used units where = 9h = P\\ = 
2^ = 1 in Eqs. (p5|)-(|32|). Near the superfluid transition 
temperature these results reduce to those by Lin-Liu et 
al|. 

The uniform texture is stable if the coefficient A is 
positive. The line where A vanishes in the v n —q plane is 
shown by solid line in Fig. 3. In a long interval L ^> £d 
the value of the wave vector q is not limited. This means 
that the uniform texture is stable only below the critical 
velocity v c \ defined by the conditions A = dA/dq = 0. 
The velocity v c \ is plotted as a function of magnetic field 
and temperature in Fig. |]. Note that v c \ vanishes in zero 
field for temperatures T < 0.85T C . 

The stability of small-angle helical textures is deter- 
mined by the coefficient B. The helix is stable if B > 
and unstable if B < 0. The stability as a function of T 
and H is indicated in Fig. ||. 

Helical textures with general opening angles were stud- 
ied numerically. The periodicity of the helix allows the 
numerical calculations to be limited to a single wave- 
length A = 2n/q using periodic boundary conditions for 
1 and d. In fact, making use of all the symmetries even 
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FIG. 4. The critical velocity v cl of the uniform texture 
(solid lines) as a function of the magnetic field H. The dif- 
ferent curves correspond to different temperatures with in- 
tervals of 0.05T C . The regions of stability and instability of 
a small-angle helix are separated by dash-dotted lines. The 
dashed line denotes v C 2 at T — T c . The corresponding critical 
current j C 2 in units of puVd is indicated by dotted line. 



a quarter of A would be sufficient. We then minimize 
the free energy (^) with respect to the five angle fields in 
©"(0)- This gives the energy as a function of v n and q: 
F(v n ,q). 

For each value of v n we determine the optimum wave 
vector g op t at which the energy of the helix is minimized. 
This process is simplified by the fact that the derivative 
of F with respect to q can be obtained using the formula 



oq q 



(33) 



with F gr and J SjX defined to be the corresponding densi- 
ties / gr and j SlX integrated over the one-dimensional in- 
terval. At the optimum value q = q pt(v n ) the derivative 
given by vanishes. The dependence of q op t on v n is 
illustrated in Fig. 3. 

The stability of helical textures is determined by the 
eigenvalues of the Hessian matrix of F(v n , q): 



dv n 2 dv n dq 
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d 2 F 
\ dv n dq 



d 2 F 
dq 2 
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A texture is stable if both the eigenvalues of the Hessian 
matrix are positive, and unstable otherwise. The stabil- 
ity region is indicated by a dashed line in Fig. 3. The 
velocity where the state at optimal wave vector g opt be- 
comes unstable is defined as the second critical velocity, 
v C 2- In helical texture the current increases much slower 
with v n than in the uniform state (Fig. pi). Therefore, the 



critical current j C 2 is substantially smaller than p\\v C 2- At 
H = and T « T c we find that v C 2 ~ 1.28i>d whereas 
j C 2 ~ 1.13p||^d H (Fig. ||). In general we find essen- 
tially no field dependence of v C 2, in contrast to v c i which 
vanishes at if = when T < 0.85T C . The temperature 
dependencies of v&, v C 2 and j C 2 in the high field limit 
are presented in Fig. ||. At high temperatures T > O.ST c 
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FIG. 5. The critical velocity as a function of temperature 
T for H 3> Hd- The different curves correspond to instabil- 
ities of the uniform texture v c i, the helical texture v C 2, the 
locked soliton u c ls, and the unlocked soliton v c \js- Helical 
textures are stable in the region between v c i and v C 2- The 
dotted line is j C 2 in units of p\\Vd. 



there is no stable helical texture and thus v c \ and v C 2 
coincide. At lower temperatures v C 2 is seen to grow dis- 
tinctly above v c i. Below 0.5T C helical textures again be- 
come unstable (Fig. [|). 

The numerical calculation of the Hessian matrix is sim- 
plified by the fact that both first derivatives ([l3]) and (|3^ ) 
are easily available. Thus the calculation of the Hessian 
at point (v n , q) requires texture minimizations only at 
three points: (v n ,q), (v n ,q + Aq) and (v n + Av n ,q). In 
order to minimize errors the number N of discretization 
points within a wave length A was kept constant and the 
discretization length Ax = X/N was varied instead. 

For all helices the opening angle grows continuously 
from zero with increasing v n , see Fig. ^| The largest sta- 
ble values for the opening angle found in the simulations 
were ~ 60 degrees. We have made numerical simulations 
in an interval containing several wave lengths of the he- 
lix, L ^> A. When the limit of stability is exceeded, it 
seems that the number of windings of the helix changes 
if a stable texture is possible at a given v n . Otherwise, 
the instability seems to lead to the growth of the opening 
angle. Most likely this leads to formation of continuous 
vortex structures |H^ ]. However, we were not able to fol- 
low this process beyond 90° opening angles because of 
the singularity in the coordinate system (P)-(p"l]). 
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FIG. 6. The opening angles of 1 in the plane _L H (J3i , 
solid line) and || H (/3j , dashed line) and the opening angle 
of d in the plane _L H (/?<j , dash-dotted line) as functions of 
driving velocity v n with q — q op t, T — 0.6T C and H 3> 



V. SOLITON TEXTURES 

In the previous section we studied the case where a 
flow was applied to an initially uniform texture. Here 
we investigate some cases where the initial state is inho- 
mogeneous. Let us consider a texture where 1 changes 
from the direction — x to x, as depicted in Fig. [?]. Such a 
texture has the property, which follows directly from the 
definition (0), that if it is rotated around x by angle 9, 
the phase difference A<£> changes by 29. Thus if nothing 
prevents the rotation of the texture, the critical current 
vanishes and the supercurrent is always dissipativc. The 
presence of magnetic field perpendicular to x prefers to 
have d, and via the dipole-dipole energy (||) also 1, in 
the plane perpendicular to H. This gives rise to a finite 
critical velocity that we aim to calculate. 

We study two different inhomogeneous structures. The 
first is known as (dipole-unlocked) soliton |i0| . This is a 
domain- wall like object where on one side d = 1 and on 
the other side d = —1. Because the change between these 
two orientations costs dipole-dipole energy, the thickness 
of the wall is on the order of the dipole length £d- In 
the absence of flow the asymptotic directions of 1 deviate 
from the normal ±x of the wall by an angle that de- 
pends on the temperature. When a small flow is applied 
perpendicular to the wall, the anisotropy of the kinetic 
energy cx — (l-v s ) 2 forces the asymptotic directions of 1 to 
±x. The calculated structure of the soliton is presented 
in Fig. |. 

The second structure we study could be called a dipole- 
locked soliton. The dipole-locking means that d(x) « 
everywhere. The region where 1 varies has finite 



(a) v n = 




(b) v n >0 




FIG. 7. Distribution of the vector 1 in a one- dimensional 
soliton structure for v n = (a) and for v n ^ in a transverse 
magnetic field H (b). The flow tends to detach the texture 
from the plane _L H. 



length only in the presence of flow. The calculated struc- 
ture of the locked soliton is presented in Fig. Q 

For both soliton structures the flow makes (3 to deviate 
from ir/2. When the flow is further increased, the struc- 
tures become unstable against unlimited winding around 
the flow direction. The critical values of v n are denoted 
by v c ls for locked soliton and w c us for unlocked soliton. 
The critical velocities are plotted in Figs. | and 0. The 
unlocked case has previously been studied by Vollhardt 
and Maki Q at T w T c using a variational approach. 
Our calculations give a much lower critical velocity than 
theirs. We note that the same process that leads to v c us 
also determines the critical velocity of a vortex sheet 
|pi|,p2[ . The locked soliton has previously been studied 
only for H = or H || v s , where the critical velocity 
vanishes and only dissipative state exists |p|,^0|-^2| . 

In the simulations the length L of the computational 
region has to be chosen large in comparison to ^d^d/i'n- 
The fast variation in the unlocked soliton sets an upper 
limit for the discretization length, which was typically 
chosen as 0.1£d- If L is not very long, the correct pro- 
cedure is to extract the critical current j c ls (and j c us) 
from the numerical calculation and then find the critical 
velocity using v cLS = jcLs/p\\ (and v cVS = j c vs/p\\)- 
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FIG. 8. The structure of the unlocked soliton at T = T c , 
H 3> Hd and v n = 0.15i>d- The angles a, j3, and <j> are defined 
in Eqs. (Jq|) — (J]~o|) and ip — tt/2 because of the high-field limit. 
The fast change within \x\ < £d is caused by the dipole un- 
locking in the soliton, and the slower change outside is caused 
by the anisotropy of the flow energy favoring 1 = ±x. 



FIG. 9. The structure of the locked soliton at T = T c , 
H S> Hd and v n — 0.8fd. The angles a, /3, and (j> are defined 
in Eqs. (Jlo|) and ip — tt/2. The length scale is purely 
determined by the flow velocity, which here is much larger 
than in Fig. H 



VI. CONCLUSIONS 

We have studied different ID textures and determined 
their stability against ID perturbations. It is reasonable 
to ask if limiting to ID perturbations is sufficient to de- 
termine the local stability. Namely, we know at least one 
situation in 3 He-A where this is not the case: uniform 
1 _L H || x. Here the ID model above gives stability until 
V c i = v c2 = y / 9d/{p± - P\\) (for H > Hi), but allowing 
v s to deviate from the direction of the applied phase dif- 
ference gives instability at v n that is by factor y/p\\/ p±_ 
lower Q. This situation differs, however, from the ones 
studied in the previous sections. For a homogeneous tex- 
ture with H _L x there exists a strict proof that only 
ID perturbations are relevant || . For helical and soliton 
textures the ID model allows a natural decay mechanism 
for the current, and we are not aware of any mechanism 
that could give a lower critical velocity. Therefore we 
believe that other than ID perturbations are unimpor- 
tant for the helical and soliton textures studied above, 
although a strict proof remains open. 

Measurements of the critical velocity are done by Ruu- 
tu et al jlB) . They study 3 He-A in a circular cylinder that 
is rotated around its axis. In the vortex-free container the 
relevant driving velocity v n = QR, where fl is the angular 
velocity and R the radius of the cylinder. Because R 3> 
£d, the flow near the cylindrical wall is one-dimensional to 
a good approximation. The magnetic field along the axis 
of the cylinder corresponds to transverse field relative to 
the flow along the whole perimeter of the cylinder, and 
thus the calculations presented above should apply to this 
case. If the field is perpendicular to the axis, all possible 



angles exist between the field and flow. In this case the 
critical velocity of the "uniform" texture is determined 
by the orientation H || v s ||, which gives a lower value 
of v c % than H _L v s . 

Ruutu et al find a considerable spread in the critical 
velocities. On one hand, our largest calculated values of 
v c correspond to the instability of the helical texture and 
coincide relatively well with the largest values observed 
in the experiments. On the other hand, the lowest mea- 
sured critical velocities can be explained by assuming the 
presence of a dipole-unlocked soliton. Quantitative com- 
parison is given by Ruutu et al ]T^ |. The comparison 
supports the basic assumption that the critical velocity in 
superfluid 3 He- A indicates an instability of the bulk, and 
it depends on the underlying texture. The bulk critical 
velocity is quantitatively better understood in superfluid 
3 He-A than in any other superfluid. 
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